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ABSTRACT 

The neutrino-driven explosion mechanism for core-collapse supernovae in its modern flavor relies 
on the additional support of hydrodynamical instabilities in achieving shock revival. Two possible 
candidates, convection and the so-called standing accretion shock instability (SASI), have been pro- 
posed for this role. In this paper, we discuss new successful simulations of supernova explosions that 
shed light on the relative importance of these two instabilities. While convection has so far been 
observed to grow first in self-consistent hydrodynamical models with multi-group neutrino transport, 
we here present the first such simulation in which the SASI grows faster while the development of 
convection is initially inhibited. We illustrate the features of this SASI-dominated regime using an 
explosion model of a 27 Mq progenitor, which is contrasted with a convectively- dominated model of 
an 8.IM0 progenitor with subsolar metallicity, whose early post-bounce behavior is more in line with 
previous II.2M0 and ISM© explosion models. We analyze the conditions discriminating between the 
two different regimes, showing that a high mass-accretion rate and a short advection time-scale are 
conducive for strong SASI activity. We also briefly discuss some important factors for capturing the 
SASI-driven regime, such as general relativity, the progenitor structure, a nuclear equation of state 
leading to a compact proto- neutron star, and the neutrino treatment. Finally, we evaluate possible 
implications of our findings for 2D and 3D supernova simulations. 

Subject headings: supernovae: general — hydrodynamics — instabilities — neutrinos — relativity 



1. INTRODUCTION 

Aspherical hydrodynamical instabilities play a 
paramount role in core-collapse supernovae. In the 
the neutrino-driven explosion mechanism in its modern 
guise, they are thought to be indispensable agents for 
enhancing the efficiency of neutrino energy deposition 
in the gain region sufficiently to eventually allow shock 
revival, but their importance does not stop there. The 
asphericities that already develop before the explosion is 
launched pre-determine the morphology of the explosion 
in the later phases, and are thus a c rucial factor for the 
mixing: instabilities in supernovae (lArnett et al.l 119891 : 
iKifonidis et al."'2003', 2006; "Hammer et al.""2010"), the de- 
velopment of pulsar kicks (Janka & Miiller 1994; Her ant 
19951: IBurrows HavesI I1996t iScheck et al.l 12004 12006 
Wongwathanarat et al.l 120101 : iNordhaus et ^11120108 



2012 ). and the remnant structure (jOawrvszczak et aLl 



20ld) 



In the purely hydrodynamical case (i.e. in the ab- 
sence of magnetic fields), two instabilities that influ- 
ence the dynamics in the supernova core in the pre- 
explosion phase have been identified. Already in the 
1990s it was recogniz ed that the gain layer can be un- 
stable to convection ( Bethel [l990f ) as the neutrino heat- 
ing of the material that is advected down to the cool- 
ing region close to the proto-neutron star surface pro- 
duces a negative entropy gradient. This was confirmed 
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by the fi rst generation of ni ulti-dimensiona l sup e rnova 
models (iHerant et al.l 119921: IB urrows & F r^^ [T99I 
'Herant et al."1994: 'Burrow s et al.hl995£ Janka fc Miillerl 
1996; Miiller & Janka 1997). The latter indeed showed 
the development of violent convective overturn in the 
heating region and suggested that the convective ex- 
change of strongly neutrino-heated material from close 
to the proto-neutron star with cool post-shock mat- 
ter in combination with the longer exposure of the ac- 
creted material to neutrino heating (jBuras et al.ll2006al : 
Murphy & Burrows 2008) can provide a powerful means 
of increasing the efficiency of neutrino heating, boost- 
ing the post-shock pressure, and eventually reviving the 
stalled shock. 

Another inst ability of a qu it e diff erent nature was 
discovered by iBlondin et al.l (j2003f ). who observed 
that a standing accretion shock as encountered in 
core-collapse supernovae may be unstable to large-scale 
i = 1 and i = 2 oscillation modes even in the absence 
of neutrino heating, a result that has been confirmed 
by further hydrodynamical sim ulations in a number of 
follow-up studies in 2D and 3D (Blondin & MezzacaDDal 
20061 : [Ohnishi et al. 2006; Blondin & MezzacaDDsy 



2007[ IScheck et al.l l20i 



Fernandez fc ThompsonI l2009b: Fer nande: 



120091; 
Lin- 
ear stability analyses provided a sound theoretical 
framework for the understanding of this "standing ac- 
cretion shock instability" and identified an amplification 
cycle of entropy /vorticity and acoustic perturbations 
between the shock and the proto-neutron star s urface as 
the underlying; mechanism , for the instabilitv (iFoglizzo 
2002[ iFoghzzo et all I2006L [20071: lYamasaki fc Yamada 



2007n. While IBlondin fc Mezzacappal (j2006D and 



LamingI ()2007[ ) have argued for the possibility of a 
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purely acoustic cycle, iGuilet fc Foglizzd (|2Q12f ) have re- 
cently given strong arguments for an advective- acoustic 
cycle. Similar to convection, strong SASI activity may 
improve the efficiency of neutrino heating by increasing 
the average shock radius and theref ore the res i dence 
time of matter in th e gain region (Scheck et all [2008[ 
iMarek Jankal 120091 ). 

While convection and the SASI are well understood 
in the linear regime, where they can be clearly differ- 
entiated, e.g. on the basis of the dominant wavenum- 
ber, w hich is typically higher for the convective modes, 
see ( Foglizzo et al.ll2006[ ) and by means of their growth 
behavior (monotonic vs. oscillatory). The situation is 
far more complicated in the non-linear regime where 
both instabilities can interact with convection trigger- 
ing secondary shock oscillations or vice versa, and the 
presence of a strong ^ = 1 or ^ = 2 mode may 
not necessarily be indicative of the SASI. The evidence 
for the SASI as "main" instability from recent 2D ex- 
plosion mod els using ela borate m ulti-gro up transport 
(iBuras et al.l i2006a: Mar ek Jank a 2009; Bruenn et all 
[20091; 'Suwaet al. 20 10: . Miiller et al.ll2012i) as weh as from 
the ffist tentative explorations in 3D (jTakiwaki et al. 
[2012) is therefore ambiguous: While they all show the 
presence of the strong i = 1 and i = 2 shock oscil- 
lations characteristic of the SASI at late times, con- 
vection appears to be the "primary" instability in the 
sense that it always begins to grow ffist a few tens 
of milliseconds after bounce. Convection also appears 
to be the "primary" instability in the 3D models of 
iFrver fc WarrenI ()2002f ) and in recent 2D and 3D sim- 
ulations based on different implemen tations of a sim- 
ple light-bulb app roximation (Murphv fc Burr ows 2008; 
Fernandez fc Thompson 2009a; Nordhaus et al. 2010b; 
Hanke et al. 2011; Burrows et al. 2012; Murphv et al. 
20121) . The work of iHanke et al.l (j201lf ) in particu- 
lar suggests that for a certain choice of parameters in 
such approximate descriptions may little activity of low- 
i modes in 3D, while the inverse turbulent energy cas- 
cade may still help to excite \ow-£ modes in 2D even 
though convection initially grows faster than the SASI. 
Except for some gray transport models w ith fast (pre- 
scribe d) proto- neutron star contraction of iScheck et al.l 
()2008f ). the primary growth of convection seems to be 
a generic property of simulations that include neutrino 
heating in some simplified or elaborate fashion. The ap- 
parent weight of evidence has even led to the sugges- 
tion "that the SASI is at most a min or feature of super- 
nova dynamics" (^Burrows et al.ll2012f ) and that neutrino- 
driven convection do minates in core-collapse supernovae 
(jMurphv et al.ll20T2h . 

However, we shall demonstrate in this paper that 
such a conclusion is more than premature, and that un- 
der favorable conditions the SASI plays a major role 
for the dynamics in the supernova core. We present 
new simulations of (incipient) explosions for two differ- 
ent progenitor stars obtained with the general relativis- 
tic multi-group neutrino hydrodynamics code Vertex- 
CoCoNuT. Based on these, we show the existence of a 
"SASI-dominated" regime distinct from a "convection- 
dominated" regime for the growth of hydrodynamical 
instabilities as a separate route towards the explosion. 
After briefly introducing the numerical methods and the 
progenitors in Section [21 we discuss the morphology of 
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Figure 1. Density profiles for model u8.1 (black solid line) and 
model s27.0 (red, dashed) covering the inner 10, 000 km. For 
comparison, profiles for the 11.2M0 progenitor sll.2 (blue, dash- 
dotted) of Wooslev et al. (2002) and the ISM© model sl5s7b2 
(green, dotted) of Wooslev Weaver (1995) are also shown. Col- 
ored circles denote the edge of the iron core and the boundary 
between the silicon shell and the oxygen-enriched silicon shell. 

the multi-dimensional flow prior and during the incipi- 
ent explosion and quantitatively analyze the growth of 
the SASI and convection in the different regimes in Sec- 
tion [3l The implications of our findings are discussed in 
Section [H 

2. NUMERICAL METHODS AND MODEL SETUP 

We perform axisymmetric (2D) core-collapse super- 
nova simulations with the general relativistic neutrino 
hydro dynamics code Vertex-CoCoNuT (Miiller et al.[ 
[2Q1Q[ ). The hydrodynamics solver CoCoNuT is a 
time-explicit, directionally unsplit implementation of the 
piecewise parabolic method (PPM) for spherical polar co- 
ordinates using a hybrid HLLC/HLLE Riemann solver, 
and relies on the extended conformal flatn ess approxima- 
tion (xCFC. ICordero-Carrion et al.[[2QQ9'[ ) for the metric 
equations. CoCoNu T is coupled to the neutrino trans- 
port module Vertex (|Rampp fc Jan ka 2002), which em- 
ploys a variable Eddington factor technique to solve the 
moment equations for the neutrino energy and momen- 
tum density with a closure provided by the formal solu- 
tion of a simplified Boltzmann equati on. In 2D, we resort 
to the ^^ray-by-ray -plus" approach (jBuras et al.[ [2006"bl : 
[Bruenn et al.[[2006[ ) , which assumes that the neutrino dis- 
tribution function is axially symmetric around the radial 
direction (and hence implies a radial flux vector). Ver- 
tex includes the velocity- and metric-dependent terms 
as well as non-isoenergetic scattering processes, thus cap- 
turing the full complexity of multi-group transport. For 
the interactions between neutrinos a nd the matter we us e 
an up-to-date set of opacities (see iMiiller et al.l [20121 ). 
For details about the implementation and a discussion 
of the accuracy of our app roach, we r efer t he reader to 
IMiiller et al.l (|20To[ ) and .Miiller et al.l (|20T2l) . 

In this study, we consid er two progenitors, na mely the 
27Mq star model s27.0 of iWooslev et al.l (|2002D with so- 
lar metallicity, and an 8.IM0 model u8.1. Th latter 
is just above the critical mass for the formation of an 
iron core a nd has an initial m etallicity of approximately 
Z = 10-"^ (|Heger et al.ll2012D P[ The structure of model 
u8.1 should be typical of such low-mass supernova pro- 

4 A model of S.OM© and Z = 10"^ did make an AGB star 
instead. 
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Figure 2. Maximum, minimum (solid lines) and average (dashed 
lines) shock radius for model s27.0 (black lines) and u8.1 (red). 

genitors, and is dramatically different from stars even a 
few tenths of a solar mass more massive. The density pro- 
file more closely resembles an asymptotic giant branch 
(AGB) star, or an electron-capture supernova progeni- 
tor, than a typical massive star structure in that it has a 
low-density "envelope" (p < 1 g cm~^) directly on top of 
a dense core {p > 10^ g cm~^) of l.SSM©, with a tran- 
sition region of only O.OSM©, mostly the carbon layer, 
in between. Outside the 1.26Mq iron core are layers of 
silicon (out to l.SOM©), oxygen (to l.SGM©), neon, and 
carbon, with implosive energy generation due to silicon, 
oxygen, and neon burning as high as 10^^ erg g~^ s~^. 
Note that the structure of model u8.1 is specific to such 
low-mass supernova progenitors, and not owing to the 
initial metallicity of the model; a different initial metal- 
licity would only change the location and extent of the 
mass range between the AGB channel and the "normal" 
channel of iron core evolution familiar from more massive 
stars. 

Model s27.0, by contrast, has a more massive and less 
compact iron core of l.SM© embedded in a thick silicon 
shell that reaches out to I.68M0, where the transition to 
the oxygen-enriched silicon shell is located. Compared 
to model u8.1, the density drops far less rapidly outside 
the iron core. In order to better illustrate the different 
density structure of the two models, we show density 
profiles of the progenitors in Figure [TJ 

We use a numerical grid of Ur x no = 400 x 128 zones 
with non-equidistant radial spacing for both progeni- 
tors. Model s27.0 was simulated using t he equation of 
state (EoS) of Lattimer & Swestyl (|1991l ) with a value 
for the bulk incompressibility modulus of nuclear matter 
of K = 220 MeV (LS220), while the softer variant with 
K = 180 MeV (LS180) was chosen for model u8.1. For 
a discussion of the validity of the latter EoS for small- 
mass (baryonic mass < 1.5Mq) proto- neutron stars de- 
spite its marginal inconsistency with the 1.97M(7) neutron 
star found by Demorest et al.j (j2010l) . see Miiller et al.l 
(^2012). Specifically, the mass-radius relation for hot and 
cold neutron stars is very similar for neutron stars well 
below the mass limit. As a consquence, both equations 
of stat e yield hardly any difference during the a ccretion 
phase (jSwestv et al.lll994l : [Thompson et al.ll2003ft . 

3. RESULTS 

Superficially, model s27.0 and u8.1 might appear to 
evolve in a very similar fashion: Roughly around 120 ms 
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Figure 3. The runaway criterion radv/'?"heat for model s27.0 
(black) and model u8.1 (red). Both Tadv cind Theat ^-^e evaluated as 
in Miill er et al.i (^2012). Note that the curves have been smoothed 
using a running average over 5 ms. 

after bounce the average shock radius starts to increase, 
and by 200 ms the shock is already expanding rapidly, 
although model s27.0 evidently lags behind u8.1 a little 
(Figure [2]). Especially during the later phases, the shock 
becomes strongly deformed with a ratio rmax/^min of the 
maximum and minimum shock radius on the order of 
2 ... 3. Both models seems to provide similar examples 
for an explosion at a relatively early stage. 

However, this appearance is deceptive: A hint at 
more profound differences between s27.0 and u8.1 is al- 
ready furnished by the critical ratio Tadv /''"heat of the 
"advection" or "residence" time-scale and the heating 
time-scale for the material in the gain region, which 
serves as an indicator for an explosive runaway due 
to neutrino energy deposition (for Tadv / The at > f ^ see 
Janka 2001; ThomDsoneFaD [2005; Buras e tall l2006al : 
iMurphv fc Burrowsi i20Q8l : pernande^ i2012|). Figure [3] 
shows that model s27.0 approaches the critical thresh- 
old much later than model u8.1, i.e. at roughly ~180 ms 
instead of ~110 ms. Nevertheless, even though the run- 
away condition is not yet met, the shock already expands 
considerably before that time in model s27.0. This sug- 
gests that at least for the first ~ 180 ms there may be a 
driving agent other than neutrino heating that is respon- 
sible for pushing the shock outwards. One should bear in 
mind, however, that it is not completely clear for which 
value of Tadv /Theat ouc could already expect a noticeable 
expansion of the shock: Neutrino heating might drive 
considerable shock expansion even for Tadv /Theat < 1 de- 
pending on progenitor specifics. However, it seems in- 
evitable that large aspherical motions in the gain regions 
with Mach numbers on the order of ~ 1 will affect the 
structure of the accretion flow, including the shock posi- 
tion (cp. Section [XT]) . 

3.1. Growth of Instabilities 

The reason for the peculiar evolution of the 27M0 pro- 
genitor is to be sought in the strong and relatively unim- 
peded growth of the SASI as primary instability dur- 
ing the first ~ 200 ms as opposed to neutrino-driven 
convection in the 8.IM0 star - a feature hitherto not 
reported fro m full multi-group neutrino hydrodynam ics 
simulations (jMarek Jankal 120091 : iMiiller et al.ll2012f ). 

As shown by Figures [5] and IH the morphology of 
the post-shock flow in model u8.1 and model s27.0 be- 
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Figure 4. Snapshots of the evolution of model u8.1, depicting the radial velocity Vr (left half of panels) and the entropy per baryon s 
(right half of panels) 66 ms, 71 ms, 83 ms, 140 ms, 200 ms, and 241 ms after bounce (from top left to bottom right). Once the gain region 
becomes convectively unstable, small-scale plumes begin to grow (66 ms, 71 ms) and merge into somewhat larger structures on the scale of 
10° . . . 50° (83 ms). As convection becomes more violent, the deformation of the shock becomes more pronounced (140 ms), and a dipolar 
asymmetry finally develops after the shock starts to accelerate outward (201 ms, 241 ms). 
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Figure 5. Snapshots of the evolution of model s27.0, depicting the radial velocity Vr (left half of panels) and the entropy per baryon s 
(right half of panels) 74 ms, 94 ms, 106 ms, 113 ms, 159 ms, and 215 ms after bounce (from top left to bottom right). The six panels 
exemplify the growth of the SASI in the linear regime (74 ms), the development of parasitic instabilities on top of a clear i = 1 asymmetry 
of the global entropy distribution (94 ms, 106 ms), the transition to the non-linear regime with the formation of a pronounced downflow 
(106 ms to 159 ms) and the phase of vigorous shock expansion (215 ms). Note the clearly different structure in the first four panels 
compared to the corresponding panels in Figure 3] 
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Figure 6. Entropy along the north and south polar axis as a function of time for models s27.0 (left) and u8.1 (right). 
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Figure 7. Coefficients for the decomposition of the shock surface into Legendre polynomials. The top panels show the ffist three normalized 
Legendre coefficients a^/ao for model s27.0 (left) and model u8.1 (right). In order to better exemplify the growth of the £ = 1 mode, \ai \ /ao 
is also plotted on a logarithmic scale in the bottom panels. The dashed line roughly denotes the slope of the exponential growth of the 
SASI up to ~120 ms in model s27.0. 




time after bounce [s] time after bounce [s] 



Figure 8. Wavelet spectrograms of the normalized Legendre coefficient ai/ao of the / = 1 mode of the shock for model s27.0 (left) and 
model u8.1 (right). Large values of the wavelet transform are indicated by darker shades. The red curves roughly define a band in which 
an advective-acoustic instability should surface; the upper and lower curves show the (smoothed) time evolution of the advection time-scale 
Tadv cind of 2Tadv5 respectively. 
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comes quite dissimilar as soon as the gain region forms 
a few tens of milliseconds after bounce. Model u8.1 
(Figure HJ conforms to the "standard" behavior exhib- 
i ted by all recent first-princ iple explosion simulations 
dMarek & Janka"2009^ 'Miiller et al. 2012): The convec- 
tive instability develops first, and perturbations at in- 
termediate scales (corresponding to multip ole orders o f 
^ > 5 ... 10) grow fastest, cf. .Fodizzo et all 120061 . 1200/D . 
Larger plumes only form once convection has become 
vigorous and the shock has already expanded quite con- 
siderably (cp. the 11.2M0 case of Miiller et al. 2012). A 
strong non-oscillatory dipolar shock deformation appears 
at late time once the explosion is underway, and this de- 
formation seems to be driven by buoyant neutrino-heated 
plumes rather than by genuine SASI activity. 

Model s27.0 behaves in a completely different manner 
(Figure [5]). Here we first observe small, but clearly recog- 
nizable linear SASI oscillations (first panel of Figure [5]). 
Later, as the amplitude of the SASI increases, parasitic 
Ray leigh- Taylor and Kelvin-Helmholtz instabilities grow 
on top of the SASI flow, but these remain localized, ap- 
pearing and disappearing as the shock oscillates back and 
forth, and never growing into developed convection. The 
different character of the post-shock flow can also be seen 
in Figure [6l which shows the time-evolution of the en- 
tropy along the polar axis for both models and illustrates 
that - unlike model 8.1, s27.0 indeed exhibits the sloshing 
motions characteristic of the SASI. Until ~ 180 ms the 
expansion and contraction of the shock in the northern 
and southern hemispheres are nicely antisynchronized. 

The velocity fluctuations in the gain region steadily 
grow in magnitude in model s27.0, eventually leading to 
the formation of secondary shocks (see the fifth panel of 
Figure [5j). It is likely that these strong velocity fluctu- 
ations are responsible for pushing the shock further out 
because they provide an additional Reynolds stress con- 
tribution for the angle- averaged and time- averaged flow 
on top of the thermal pressure. Moreover, the thermal 
pressure itself is increased once dissipation in secondary 
shock starts to convert kinetic energy stored in the SASI 
motions into thermal energy. The dynamical importance 
of both these effects depends on the typical velocity 5v of 
the aspherical flow perturbations, which determines the 
ratio of the Reynolds stresses p5v'^ and the kinetic en- 
ergy to the thermal pressure Ptherm and energy density 
petherm, respectively: 

p^V^/Ptherm CX {l/2p5v'^) / (petherm) CX 5v'^ / cl = Ma^ 

(1) 

Once the SASI reaches the non-linear stage in model 
s27.0, the typical Mach number Ma (determined from 
angular root mean square averages of the lateral veloc- 
ity component) exceeds > 0.5 in large parts of the gain 
region and peaks at 1 near the average shock radius, 
clearly indicating the dynamical importance of aspheri- 
cal instabilities in that model. In that sense, model s27.0 
may provide an example for S ASI-driven shock expansion 
prior to the explosion as opposed to the neutrino-driven 
shock expansion that roughly sets in as radv/'?"heat exceeds 
unity. 

We emphasize that these mechanisms for shock ex- 
pans ion are, of course, n ot specific to the SASI; in- 
deed i Murphv et al.l (|2012) demonstrated that the inclu- 
sion of the Reynolds stresses in the shock jump condi- 



tions largely accounts for the increased shock radii in 
their convective models. However, while the immediate 
cause for the expansion of the shock is similar in model 
s27.0 to convectively- aided explosion models like u8.1 or 
those of Murphy e t al.l ((2012), the ultimate source for 
the energy of aspherical motions is different: As long as 
the SASI remains the dominant instability, it essentially 
feeds directly on the kinetic energy of the accreted matter 
whereas convection is (indirectly) powered by neutrino 
heating. 

3.2. Evidence for the SASI - Quantitative Analysis 

For model s27, we thus have a SASI-like flow morp hol- 
ogy r eminiscent of si mulations with no jBlondin et al.l 
l2003f ) or suppressed ()S check et al.l 12008) neutrino heat- 
ing. The associated aspherical velocities apparently be- 
come large enough (i.e. on the order of the speed of 
sound) to be dynamically relevant. But model s27 fur- 
nishes further evidence pointing to the SASI as the ini- 
tially dominant instability besides the qualitative flow 
morphology. 

A mode analysis of the shock surface, based on a 
decomposition into Legendre polynomials, confirms the 
dramatic differences in the growth of hydrodynamical in- 
stabilities (Figure [7|). The expansion coefficients are 
given by 

at = / rsi,{0)Pt{cosO)smOdO, (2) 



and the first three coefficients ai, a2, and as (normal- 
ized to the angle- averaged shock position (rgh) = <^o)- 
The convectively-dominated model u8.1 (right panel of 
Figure [7j) is characterized by relatively weak shock oscil- 
lations until ~ 130 ms when the explosive runaway due 
to neutrino heating is already underway. Afterwards, a 
pronounced dipole (^ = 1) and quadrupole {i = 2) mode 
appear, but there is no trace of periodic oscillations at 
any stage. Apparently, the growth of the low-^ modes is 
a stochastic process driven by the convective plumes in 
this case. Later in the explosion, the geometry becomes 
more stable with a big buoyant bubble in the southern 
hemisphere. The dominance of such large-scale bubbles 
at late times is in agreement with analytic estimates for 
extended convective regions with a small, almost po int- 
like h eating source (Chandrasekhar 1961; Foglizzo et alj 
f200l . 

By contrast, model s27.0 (left panel of Figure [7|) ex- 
hibits shock oscillations with a well-defined periodic- 
ity for the i = 1 mode as expected for a relatively 
unperturbed SASI (cf. the numerical experiments by 
IS check et al.l 12008*) . The expected exponential growth 
of the dipole amplitude \ai\/ao from ~ 0.02 to ~ 0.4 
during the linear phase can be observed for several os- 
cillation periods, and a growth rate cjsasi ^ 45 s~^ 
can be deduced. The quadrupole amplitude {£ — 2) 
increases in a non-oscillatory manner after ~ 90 ms, 
which may be a result of the transition to the non-linear 
regime. Moreover, we expect the oscillation period to 
scale with the advection time-scale by a factor between 
unity and < 2 if an advectice-acoustic cycle operates 
i n model s27.0 f Foglizzo etap 120071 : IScheck et al.|[2QQ8l : 
iGuilet fc Foglizzg .2012.) . This confirmed by a wavelet 
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Figure 9. Lateral velocity vq in the equatorial plane as a function of time and radius for model s27.0 (left) and model u8.1 (right). 
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Figure 10. Evolution of the stability parameter x for model s27.0 
(black) and model u8.1 (red). The threshold value of x ~ 3 for 
convection is reached about 60 ms after bounce for model u8.1, 
whereas convective stability applies for s27.0 until 120 ms when 
the SASI has already grown to the non-linear regime. Note that 
X becomes less meaningful as the post-shock flow develops strong 
asphericities. For this reason, e.g., the fact that x drops below 
three at f ~ 0.08 s once convection has become vigorous does not 
indicate declining convective activity. 




0.1 0.15 
time after bounce [s] 

Figure 11. Advection time-scale Tadv (solid lines) and mass accre- 
tion rate M (dashed lines) for model s27.0 (black) and model u8.1 
(red). The time-scales are evaluated as described in Miiller et al. 
(^2012 ), and the curves have been smoothed using a running average 
over 5 ms. 



analysis of the 1 = 1 co efficient of the shock posit ion us- 
ing the Morlet wavelet (|Torrence fc CompQlll998[ ) shown 
in Figure [H Even after ^150 ms, when the underlying 
linear SASI eigenmode can no longer be easily discerned 



in the post-shock flow (see the last two panels in Figure [5] 
andH]), the advection time-scale apparently still sets the 
oscillation period. This is in stark contrast to model 
u8.1, for which the wavelet spectrogram never shows a 
lot of power in the expected "SASI frequency band" . The 
antisynchronized variation of the shock radius along the 
north and south polar axis (left panel of Figure [6j) that 
can be discerned until ~ 180 ms also points to a global 
low-^ instability; such a correlated variation of the shock 
radius in opposite directions is probably hard to reconcile 
with a stochastic forcing of low-^ modes by convection. 

Up to ~ 150 ms, a closer look at the post-shock velocity 
field even reveals the advective part of the SASI ampli- 
fication directly. Figure [9] clearly shows the advection 
of coherent perturbations in the lateral velocity (consis- 
tent with the i = 1 sloshing mode) from the shock to 
the deceleration region at the proto-neutron star surface 
for about 100 ms. Around ~ 150 ms, the perturbations 
become somewhat crinkled due to the action of parasitic 
instabilities, but a temporal quasi-periodicity can still be 
recognized. In model u8.1, we see only very faint traces 
of such coherent perturbations prior to the onset of con- 
vection. 

There are thus good reasons for classifying model s27 
as "SASI-dominated" and model u8.1 as "convectivion- 
dominated" at least for first ~ 150 ms. The late-time 
behavior around and after shock revival is probably more 
complicated in both cases. During this phase, neutrino 
heating evidently plays a major role for the dynamics and 
ultimately powers the explosion, and buoyancy-driven 
instabilities certainly become important in the process. 
However, this does not imply that steady- state convec- 
tion alone furnishes an appropriate picture for the post- 
shock fiow at this stage. The simulations rather sug- 
gest a complicated interplay between neutrino-heated, 
buoyancy-driven bubbles, internal shocks and acoustic 
waves, supersonic downflows, and an aspherical shock 
that creates large non-radial velocities in the post-shock 
region (bottom panels of Figures H] and [5]) . The charac- 
terization of this phase will require a better understand- 
ing of the interaction of the different hydrodynamical 
instabilities in the fully non-linear phase. 



3.3. Conditions for SASI- Dominated and 
Convectively- Dominated Flow 
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Why are the two models s27.0 and u8.1 dominated 
by different hydro dynamical instabilities, and why, in 
particular, is convection suppressed for the 21 Mq pro- 
genitor despite neutrino heating in the gain layer? The 
key for understanding the differences is connected to the 
fact that the presence of a negative entropy gradient in 
the heating region is not sufficient for the onse t of con- 
vection as pointed out by iFoglizzo et al.l (|2QQ6l ) . Under 
the assumption of a stalled accretion shock, any con- 
vective perturbation is advected out of the gain layer 
within a finite time, and convection can develop only 
if the perturbation is amplified sufficiently within this 
tim e frame. Using a mor e rigorous mathematical analy- 
sis, IFoglizzo et al I (|2QQ6[ ) derived a growth parameter x 
that is defined in terms of the Brunt- Vaisala frequency 
cjbv and the spherically averaged advection velocity {vr) 
as follows, 

Irm^ 

where the integral runs from the gain radius Vg 
to the average shock radius rgh- Note that only 
the region where cjgy < indicates local instabil- 
ity in the fluid fr ame contributes to the integral. 
IFoglizzo et al.l ()2QQ 6^ find a threshold condition of x ^ 
3 for convective instability in the gain region, which 
has been confirmed by parametrized as well as first- 
princ i ple simulations (Buras e t al. 2006a; Scheck et al. 
120081 : iFernandez Thompsonl l2009bt ). We empha- 
size that the criter ion X ^ 3 has been derived by 
IFoglizzo et al.l (|2006l ) for the linear regime and may no 
longer be useful for determining the presence or absence 
of convective instability once significant non-spherical 
perturbations develop in the post-shock region. 

The time evolution of x for our two models is shown 
in Figure [TOl Clearly, the parameter x indicates that 
convection should set in at ~60 ms for u8.1, which is ex- 
actly what we observe (Figure H]). For s27.0, on the other 
hand, x is very low at this junction, and only reaches the 
critical value x = 3 some 120 ms after bounce, when the 
flow is already strongly aspherical due to the activity of 
the SASI (Figure [7]) and x i^iay no longer be a reliable 
measure for convective instability. 

The conditions for effective SASI growth are quite dif- 
ferent. Here, a short advection time-scale Tadv is helpful: 
In the picture of the advective-acoustic cycle, the linear 
growth rate a;sASi is given in terms of the cycle efficiency 
Q and the duration Tcyc of the cycle by 

^SASI = ^^, (4) 

Tcyc 

as shown by IFoglizzo et al.l (|2006[ ). A smaller shock 
stagnation radius and advection time-scale also imply a 
shorter cycle duration Tcyc, and hence conditions favor- 
able for a more rapid growth of the SASI0 Since Tadv 
is shorter by a factor of 2 ... 4 in model s27.0 compared 
to u8.1 (Figure [TT]), considerably more vigorous SASI 
activity is to be expected. Moreover, convection could 
possibly destroy the coherence of the entropy, vorticity, 

^ Qualitatively, the same dependence would be expected for a 
purely acoustic cycle. Whether the SASI grows due to advective- 
acoustic or purely acoustic feedback is therefore irrelevant for our 
discussion. 



and acoustic wave s involved in the amplification cycle 
(jGuilet et al.ll2010D even if it develops as a primary and 
not as a parasitic instability so that the very absence of 
convection for the 27 Mq progenitor may also be a crucial 
factor for violent SASI activity. 

Given these conditions, the predominance of the SASI 
and convection in model s27.0 and model u8.1, respec- 
tively, can be well accounted for, but the growth condi- 
tions can be further connected to the progenitor struc- 
ture. The shorter advection time-scale in s27.0 is a direct 
result of a significantly higher accretion rate (by a factor 
of ~ 4 at 150 ms), which is in turn a direct consequence 
of the presence of a dense and rath er massive silicon 
shell on top of the iron core (see, e.g.. IWooslev fc Hegerl 
'2012' for the relation between the progenitor structure 
and the time-dependence of M). The higher mass ac- 
cretion rate leads to a smaller shock stagnation radius 
(Figure [2]) after 30 ms post-bounce and thus to a more 
narrow post-shock layer as well as higher pre- and post- 
shock velocities, which implies a reduction of the accre- 
tion time-scale. Particularly at later stages, the stronger 
gravitational field of the proto- neutron star, whose bary- 
onic mass reaches some l.GSM© compared to 1.36M0 for 
u8.1 also contributes to the shorter advection time-scale. 
According to Equation ([3]), this will reduce the value of 
X and therefore inhibit the development of convection. 

Potentially, the larger advection velocities in the gain 
region could be compensated by a higher value of ujbv 
in Equation (|3]), which might result from the stronger 
neutrino heating that may be expected because of an en- 
hanced accretion luminosity (and hence stronger heating) 
in model s27.0. Larger neutrino luminosities and mean 
energies are indeed observed for model s27.0 (Figure [T2|). 
but what does this imply for the instabilities in the gain 
region? 

To decide this question, we estimate the dependence of 
the unstable non-adiabatic gradient, e.g., of the internal 
energy density e, on the neutrino emission and the heat- 
ing efficiency by means of simple zeroth-order approxi- 
mation for the solution of the internal energy equation 
for an infalling mass shell: 

— (X (X — — (X — 

V y heating ^sh - Tg Mg (r^h - Tg) M (fgh - Tg) 

(5) 

Here, is the mass-specific neutrino heating rate, which 
is approximated in terms of the volume- integrated heat- 
ing rate Qj^ and the mass in the gain region Mg ^ Mradv 
Equation (|5]) shows that stronger heating might be bal- 
anced by a larger M, but as Qj^ still depends on M both 
through the post-shock stratification and the accretion 
luminosity Equation dU does not help very much. 

It is more useful to reformulate Equation dU in terms 
of two efficiency parameters 77heat and 77conv for the neu- 
trino heating and for the conversion of gravitational en- 
ergy into electron (anti-) neutrinos. The familiar defi- 
nition of the heating efficiency ( Marek fc Jankal 120091 : 
iMiiller et al.l l2012f ) in terms of the volume-integrated 
heating rate Qheat in the gain region and the electron 
neutrino and antineutrino luminositie^ and L^^^, 

^ We use the total angle- integrated neutrino energy flux mea- 
sured at inflnity (FigureEJ for -^i^e ^^<1 ^^^e • Note that ?7heat could 
be interpreted as an effective (non-Rosseland) frequency-averaged 
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Figure 13. Heating efficiency ?7heat (solid lines) and conversion 
efficiency ?7conv (dashed lines) for s27.0 (black) and u8.1 (red). All 
curves have been smoothed using a running average over 5 ms. 



reads 



^heat 



(6) 



Figure 12. Neutrino luminosities (left panel) and mean energies (right panel) for s27.0 (black) and u8.1 (red). Solid, dashed, and dash- 
dotted lines are used for Ue, i^e, and ^n/ri respectively. Note that we plot angle- aver aged quantities extracted at a fiducial observer radius 
of 400 km. 

through the gain layer (Figure [TT]), which is compensated 
by a larger mass in the gain region at later times. Par- 
ticularly large differences are observed for T^conv, which 
is larger by a factor of 2 ... 4 in model u8.1 prior to the 
explosion and is therefore by far the most relevant term 
in Equation ([8]). At early times, the rapid post-shock 
advection in model s27.0 also results in somewhat ineffi- 
cient cooling, whereas neutrino diffusion from the deeper 
regions of the proto-neutron star starts to enhance the 
electron (anti-) neutrino luminosity rather early for the 
8.IM0 star and soon outweighs the accretion luminosity. 
As a consequence, the relative difference in luminosity 
between the two models (Figure [12]) is much smaller than 
the relative difference of the accretion rates (Figure [TT]) . 
The two terms in Equation ([8]) that could enhance the 
unstable gradient in s27.0 compared to u8.1 are the differ- 
ence of the shock and gain radius in the denominator and 
the compactness parameter M / R. However, the effect of 
the different value for rgh — is more than compensated 
for by different domain of integration in Equation (|3]) 
for the parameter x-> ^^^^ the compactness M/R is only 
higher by - 25%. 

Our analysis thus provides the clue for understanding 
why the larger luminosities in model s27.0 are not helpful 
for fostering the growth of convective instabilities: The 
decisive quantity is the ratio between the luminosity and 
a fiducial accretion luminosity GMM/R (rather than the 
luminosity itself), and a high accretion rate will decrease 
rather than increase this "conversion efficiency" . 

4. DISCUSSION AND CONCLUSIONS 

We presented new two-dimensional general relativis- 
tic explosion models for a metal-poor 8.IM0 and solar- 
met allicity 27 Mq progenitor. Our successful explosion 
models for these stars exemplify two possible regimes for 
the growth of hydrodynamic instabilities in core-collapse 
super novae. While the 8.IM0 star conforms to the fa- 
miliar picture of neutrino-driven convection growing first 
and giving rise to low-^ shock oscillations afterwards as 
shock exp ansion sets in (se e, e.g., the II.2M0 and 15 M^ ;^ 
models of iMarek Jankal l2QQ9 and IMiiller et all 120121 ). 
the SASI is the primary instability for the 27 Mq pro- 
genitor. In this model, we can clearly identify the SASI 
and its characteristic features during the linear growth 
phase, and find that the flow remains distinctively differ- 



As a measure for the efficiency of the conversion of grav- 
itational energy of the accreted material into neutrinos 
we define r^conv in terms of L^^, L^^, and M as 



Vc 



{L^^^LpjR 



(7) 



GMM 

Small values of ?7conv < 1 indicate an inefficient conver- 
sion of gravitational energy into accretion luminosity due 
to the fast advection of the accreted material through 
the cooling layer, whereas r^conv > 1 indicates an effec- 
tive conversion and an additional enhancement of the 
electron (anti-)neutrino luminosity due to the diffusive 
neutrino flux from the deeper layers of the proto-neutron 
star. The unstable gradient manifestly depends on those 
two efficiencies: 



oc 



heating 



con\ 
^(^sh - rg) 



(8) 



As shown in Figure [131 both reheat and 77conv remain sig- 
nificantly smaller for s27.0 than for u8.1 before the SASI 
starts to push the shock outwards at ~ 120 ms. The 
smaller heating efficiency ?7heat in s27.0 at early times is 
a direct consequence of the rapid advection of material 

optical depth of the gain region. 
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ent from the convective S.IMq model in the non-linear 
regime. This "SASI-dominated" regime is observed here 
for the first time in a fully self-consistent neutrino trans- 
port simulation; the purely hydrodynamical models of 
iBlondin etlo] (|2QQ3[ ). in which the SASI was initially 
discovered, and the SASI models of Scheck et al. (2008), 
which relied on a fast contraction of the proto-neutron 
star, may thus capture the dynam ics in the supernova 
core better than re cently suggested ([Burrows eraII[20Tl 
iMurphv et al.l 1201 2). Whether convection or the SASI 
emerges as the dominant instability evidently depends 
on the conditions in the accretion flow onto the proto- 
neutron star. Moreover, the SASI plays a major role 
during the evolution of the 27 Mq progenitor towards an 
explosion as it already pushes the shock out to an av- 
erage radius of 300 km before the runaway condition 
for neutrino-driven shock expansion is finally reached 
and parasitic convective activity becomes very strong. 
It is notew o rthy that, contrary to the e xpect ations of 
iFrver et al.l (|2012[ ) and iBelczvnski et al.l (|2011[ ). such a 
SASI-aided explosion can be initiated at a similarly early 
stage as in the convectively-dominated model u8.1, i.e. 
less than 200 ms after bounce. 

There may be a number of conspiring factors 
that allow the 27 Mq progenitor to reach the SASI- 
driven regime: Any effect that leads to higher 
post-shock advection velocities could potentially de- 
cide about the character of the dominant aspheri- 
cal instability in the pre-explosion phase. For our 
27 Mq model, the following differences to other re- 
cent two-dimension al first-principle siniulations of su- 
pernova explosions (|Bu ras et al.ll2006al: "Marek & Janka' 
| 200aiBruenn et al.l [2^9: Suwa et al. 2010; Miiller et al. 
I2OI2I ) as well as convective models based on a sim- 
ple light-bulb approximat ion (M urphv & Burrows 2008t 
[Nordhaus et al. 2010b; Ha nke et alJl201ll : iBurrows et all 
[2012; Mu rphv et al. 2012) contribute - certainly to vary- 
ing degrees - to the short advection time-scale: 

1. The iron core of the progenitor - and hence the 
proto-neutron star - is fairly massive, and is sur- 
rounded by a thick and rather dense silicon shell. 
The mass accretion rate therefore remains high 
long after the gain region has formed. This implies 
a small shock stagnation radius and a short advec- 
tion time-scale, which in turn suppresses convective 
activity (see Equation [3]) . On the other hand, the 
short advection time-scale is conducive to efficient 
SASI growth (see Equation H]) . 

2. General relativity further contributes to the short 
advection time-scal e as it leads to a more com- 
pact neutron star (Bruenn et al.|[2QQlt iMiiller et al.l 
I2OI2I : iLentz et al.ll2012f ). an "effect that is ah the 
more important as the proto-neutron star is al- 
ready quite massive to begin with. Newtonian 
models could thus systematically miss the SASI- 
driven regime. 

3. The nuclear EoS LS220 produces fairly compact 
neutron stars in ag;reement with recent radius es - 
timates (jSteiner et al.l 120101 : iHebeler et aIll2010D . 
again contribut ing to fast adve ction through the 
gain layer (cp. iMarek fc Jankal [20091 for the de- 
pendence of the accretion shock radius on the 



EoS). S ASI growth may be inhibited in si mu- 
lations (IBurrows et al . 2006; Murphv & Burr owsl 
120081 : IBurrows et al."2012: Murph v et al...2012.) us- 
ing the Shen EoS (Shen et al. 1998|) with its high 
value for the symmetry energy of 36.9 MeV. 

4. Including the full neutrino transport for all flavors 
allows the proto-neutron star to contract properly, 
whereas light-bulb models and simulations neglect- 
ing fi and r neutrinos may suppress or underesti- 
mate the contraction (and hence the advection ve- 
locities) considerably. Light-bulb models that dis- 
regard neutrino diffusion from the neutron star in- 
terior neglect this aspect of proto-neutron star evo- 
lution completely. 

With all these elements working in tandem to decrease 
the accretion shock radius prior to the onset of multi- 
dimensional instabilities, the competition between buoy- 
ancy effects and post-shock advection is heavily tilted to- 
wards the latter, res ulting in a suppression of convection 
(jFodizzo e t Id] [200 6). while the growth of the SASI is 
accelerated. The individual contribution of each factor, 
as well as certain balancing effects due to enhanced elec- 
tron neutrino and antineutrino luminosities and spectra, 
would obviously merit further investigation, and cannot 
be quantifled based on a single model. 

The realization that the SASI rather than convec- 
tion can be the primary instability under certain con- 
ditions encountered in self-consistent 2D supernova sim- 
ulations opens up several interesting perspectives: Can a 
model that is initially dominated by convection perhaps 
undergo a late-time transition to the SASI-dominated 
regime if the shock retracts sufficiently far? This could, 
for example, provide an ex planation for the lat e ex- 
plosions of a ISM^ star of IMarek fc Janka[ (|2009[ ) and 
IMiiller et al.l (|2012l) . which are associated with strong 
shock oscillations. Furthermore, the growth and satu- 
ration behavior of the SASI in 3D become a major is- 
sues: Could the parasitic instabilit ies that have been 
proposed as saturation mechanism (Guilet et alj |201Q[ ) 
be sufficiently suppressed to allow vigorous SASI activ- 
ity in 3D as well? Could even a partial suppression 
due to fast advection explain the development of low- 
i modes, which may be inhibited under other circum- 
stances becaus e of the forward turbulent cascade in 3D 
(Hanke et al. 2QTT[ )? Wih the ^ = 1 sloshing mode be 
replaced b y a spiral mode that inight spin up the neu- 
tron star ( Blondin fc Mezzacappall2007l : lFernandezll2010D 
much mor e effectively than suggested by recent sim- 
ulatio ns (Wongwathanarat et ajj 120101 : iRan tsiou e t al.l 
|2011[ )? Will the SASI saturate at lower amplitudes in 3D 
as the kinetic energy is sh ared among a larger number 
of modes as suggested bv Iwakami et al.l ((2009), or will 
the kinetic energy contained in non-radial motions even 
become larger as more modes can be excited? Will the 
growth and saturation depend on the initial seed pertur- 
bations and the rotation rate of the progenitor? All these 
questions will need to be answered not only for restricted 
parametrized setups^ but for the realistic conditions ac- 
tually encountered in self-consistent neutrino hydrody- 
namics simulations of a large variety of progenitors. The 
future exploration of the SASI and convection will also 
require better methods for discriminating between these 
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two - possibly sometimes co-existing - instabilities in the 
non-linear regime. The role of SASI and convection in 
core-collapse supernovae as well as their mutual interac- 
tion are thus bound to remain fruitful topics in supernova 
physics. 
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